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ABSTRACT 

We show that the mass-metalhcity relation observed in the local universe is due to a 
more general relation between stellar mass M^, gas-phase metallicity and star forma- 
tion rate (SFR). Local galaxies define a tight surface in this 3D space, the Fundamental 
Metallicity Relation (FMR), with a small residual dispersion of ~0.05 dex in metallic- 
ity, i.e, ^12%. At low stellar mass, metallicity decreases sharply with increasing SFR, 
while at high stellar mass, metallicity does not depend on SFR. 

High redshift galaxies, up to z~2.5 are found to follow the same FMR defined by 
local SDSS galaxies, with no indication of evolution. In this respect, the FMR defines 
the properties of metal enrichment of galaxies in the last 80% of cosmic time. The 
evolution of the mass-metallicity relation observed up to z=2.5 is due to the fact that 
galaxies with progressively higher SFRs, and therefore lower metallicities, are selected 
at increasing redshifts, sampling different parts of the same FMR. 

By introducing the new quantity fia = log(Mj,) — alog(SFR), with a=0.32, we 
define a projection of the FMR that minimizes the metallicity scatter of local galaxies. 
The same quantity also cancels out any redshift evolution up to z^2.5, i.e, all galaxies 
follow the same relation between /io.32 and metallicity and have the same range of 
values of /io.32- At z>2.5, evolution of about 0.6 dex off the FMR is observed, with 
high-redshift galaxies showing lower metallicities. 

The existence of the FMR can be explained by the interplay of infall of pristine 
gas and outflow of enriched material. The former effect is responsible for the depen- 
dence of metallicity with SFR and is the dominant effect at high-redshift, while the 
latter introduces the dependence on stellar mass and dominates at low redshift. The 
combination of these two effects, together with the Schmidt-Kennicutt law, explains 
the shape of the FMR and the role of /io.32- The small metallicity scatter around the 
FMR supports the smooth infall scenario of gas accretion in the local universe. 

Key words: Galaxies: abundances; Galaxies: formation; Galaxies: high- 
redshift;Galaxies: starburst 



1 INTRODUCTION 

Gas metallicity is regulated by a complex interplay between 
star formation, infall of metal-poor gas and outflow of en- 
riched material. A relation between magnitude and metal- 
licity was discovered in the '70 (McClure & van den Bergh 
1968; Lequeux et al. 1979), in which more luminous galaxies 
also have higher metallicities. Later on, it was understood 
that this luminosity-metallicity relation is a manifestation of 
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a more fundamental stellar mass-metallicity relation where 
galaxies with larger stellar mass M* have higher metallicities 
(Garnett 2002; Perez-Gonzalez et al. 2003; Pilyugin et al. 
2004; Tremonti et al. 2004; Savaglio et al. 2005; Lee et al. 
2006; Cowie & Barger 2008; Panter et al. 2008; Kewley & 
Ellison 2008; Hayashi et al. 2009; Michel-Dansac et al. 2008; 
Liu et al. 2008; Rodrigues et al. 2008; Perez-Montero et al. 
2009). 

The origin of this relation is debated, and many differ- 
ent explanations have been proposed, including ejection of 
metal-enriched gas (e.g., Edmunds 1990; Lehnert & Heck- 
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man 1996; Garnett 2002; Tremonti et al. 2004; Kobayashi 
et al. 2007; Scannapieco et al. 2008; Spitoni et al. 2010), 
"downsizing", i.e., a systematic dependence of the efficiency 
of star formation with galaxy mass (e.g., Brooks et al. 2007; 
Mouchine et al. 2008; Calura et al. 2009), variation of the 
IMF with galaxy mass (Koppen et al. 2007), and infall of 
metal-poor gas (Finlator & Dave 2008; Dave et al. 2010). 

Recently, evidence has been found that at high redshift 
the infall of pristine gas can have a dominant role (Bour- 
naud & Elmegreen 2009; Dekel et al. 2009; Brooks et al. 
2009; Agertz et al. 2009). The mass-metallicity relation has 
been studied by Erb et al. (2006) at z~2.2 and by Maiolino 
et al. (2008) and Mannucci et al. (2009) at z=3-4, finding a 
strong and monotonic evolution, with metallicity decreasing 
with redshift at a given mass. The same authors (Erb et al. 
2006; Erb 2008; Mannucci et al. 2009) have also studied the 
relation between metallicity and gas fraction, i.e., the effec- 
tive yields, obtaining clear evidence of the presence of infall 
in high redshift galaxies. 

Even more recently, Cresci (2010) have obtained the ev- 
idence of "positive" metallicity gradients in a small sample 
of three disk galaxies at z=3.5, with lower metallicities in 
more active regions of star formation. The presence of these 
gradients in galaxies with very regular dynamics can be un- 
derstood as a consequence of accretion of metal-poor gas, 
producing a new episode of star- formation in older, more 
metal-rich galaxies. 

If infall is at the origin of the star formation activity, 
and outflows are produced by exploding supernovae (SNe), 
a relation between metallicity and SFR is likely to exist. In 
other words, SFR is a parameter that should be considered 
in the scaling relations that include metallicity. The role of 
specific SFR (SSFR) in this content was already studied 
by Ellison et al. (2008), who presented a mild (^0.1 dex) 
dependence of metallicity on SSFR when binning galaxies 
according to their M*. Recently, also Lopez-Sanchez (2010) 
presented evidence for a link between SFR and metallicity, 
while Peeples et al. (2009) reported high values of SFR in 
a sample of outliers, toward low metallicities, of the mass- 
metallicity relation. 

To test the hypothesis of a correlation between SFR 
and metallicity in the present universe and at high redshift, 
we have studied several samples of galaxies at different red- 
shifts whose metallicity, M*, and SFR have been measured. 
In the next section we present the data samples we are 
using for this study. In sec. 3 we study the mass-metallicity 
relation as a function of SFR, and in the following section 
we introduce the Fundamental Metallicity Relation. In 
sec. 7 we discuss the physical origin of this relation. We 
adopt a ACDM cosmology with Ho—70 km/sec, 0^=0.3 
and r2A=0.7. Stellar mass M* and SFR are expressed in 
M0 and in Mq/jt, respectively. 



2 THE GALAXY SAMPLES 

2.1 z=0: SDSS 

Local galaxies are well measured by the SDSS project 
(Abazajian et al. 2009). We used the MPA/JHU catalog 
of emission line fluxes and stellar masses from SDSS-DR7 



available at http://www.mpa-garching.mpg.de/SDSS and 
described in Kauffmann et al. (2003b), Brinchmann et al. 
(2004) and Salim et al. (2007). The catalog includes 927552 
galaxies whose spectroscopic properties, such as emission 
line fluxes and spectroscopic indexes, have been measured. 

We selected emission-line galaxies with redshift between 
0.07 and 0.30 (47% of the total sample). The minimum 
redshift is set in order to ensure that [OII]A3727 is well 
within the useful spectral range, and that the 3" aperture of 
the spectroscopic fiber samples a significant fraction of the 
galaxies (3" correspond to ~4 kpc at z=0.07). 

A threshold to the signal-to-noise ratio (SNR) of Ha 
of SNR>25 was used to have reliable values of metallicity, 
while no SNR threshold is used for the other lines. Such a 
high SNR on Ha is needed to ensure that all the main optical 
lines are generally detected with enough SNR without intro- 
ducing metallicity biases. For example, the [NIIJA6584 flux 
is about 1/2 of the Ha flux at high metallicities, and about 
1/10 at the lowest metallicities sampled by SDSS. Noise and 
intrinsic dispersion of the [N11]A6584/Ha ratio could pro- 
duce very low fluxes of the [NIIJA6584 in low metallicity 
galaxies. If a threshold in the SNR of [NII]A6584 is used, 
some low-metallicity galaxies would be removed from the 
sample, while this effect would be negligible among high- 
metallicity galaxies. As a result a metallicity bias would be 
introduced, and such a bias can be important when high 
SNR thresholds are used. Also, this selection bias would be 
correlated with SFR as more active galaxies have brighter 
lines. The use of a SNR threshold of 25 on Ha means that, 
on average, the faintest [N11JA6584 lines are detected with 
SNR>2.5, but even galaxies with lower SNR for this line are 
included in the sample. This limit on SNR selects 43% of the 
remaining sample. 

We limited dust extinction to Ay < 2.5, in order not to 
deal with very large extinction corrections, and galaxies with 
Balmer decrements below 2.5 were removed. These selections 
remove 0.2% of the galaxies. Finally, AGN-like galaxies (22% 
of the sample) were excluded by using the BPT classification 
by Kauffmann et al. (2003a). 

Total stellar masses M* from Kauffmann et al. (2003b) 
were used, as listed in the same MPA/JHU catalog, with 
a correction factor of 1.06 to scale the masses down from 
a Kroupa (2001) to a Chabrier (2003) initial mass function 
(IMF). 

SFRs inside the spectroscopic aperture were measured 
from the Ha emission line flux corrected for dust extinction 
as estimated from the Balmer decrement. The conversion 
factor between Ha luminosity and SFR in Kennicutt (1998) 
was used, corrected to a Chabrier (2003) IMF. 

Oxygen gas-phase abundances were measured from the 
emission line ratios as described in Nagao et al. (2006) and 
Maiolino et al. (2008). Two independent measurements of 
metallicity are available for these galaxies, based either on 
the [NII]A6584/Ha ratio or on the R23 quantity, defined 
as R23=([OII]A3727-h[OIII]A4958,5007)/H/3. When both 
quantities are inside the useful range for metallicity cali- 
bration, i.e., log([NII]A6584/Ha)<-0.35 and log(R23)<0.90, 
we selected only galaxies where the two values of metallicity 
differ less that 0.25 dex (97% of the sample), and galaxy 
metallicity is then defined as the average of these two val- 
ues. The spread of the difference between these two esti- 
mates of metallicity is 0.09 dex (~23%), with a significant 
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Figure 1. Left panel: The mass-inetallicity relation of local SDSS galaxies. The grey-shaded areas contain 64% and 90% of all SDSS 
galaxies, with the thick central line showing the median relation. The colored lines show the median metallicities, as a function of M<^, 
of SDSS galaxies with different values of SFR. Right panel: median metallicity as a function of SFR for galaxies of different M, . At all 
Mi, with log(M*)<10.7, metallicity decreases with increasing SFR at constant mass . 



albeit small systematic difference of 0.05 dex (~12%) with 
the value from R23 systematically higher than that derived 
from [NII]A6584/Ha. This smaU difference is likely to be 
due to the different sample used here and in Maiolino et al. 
(2008), which use a SNR threshold of 10 on the flux of each 
line. This may introduce a small bias in the calibrations rel- 
ative to our sample. 

The final galaxy sample contains 141825 galaxies. 



2.2 z=0.5-2.5 

Many galaxies has been observed at high redshift and these 
data can be used to study the evolution of metallicity 
with respect to the other properties of galaxies. We ex- 
tracted from the literature three samples of galaxies at 
intermediate redshifts, for a total of 182 objects, having 
published values of emission line fluxes, Af*, and dust ex- 
tinction: 0.5<z<0.9 (Savagho et al. 2005, GDDS galaxies), 
1.0<z<1.6 (Shapley et al. 2005; Liu et al. 2008; Wright 
et al. 2009; Epinat et al. 2009), and 2.0<z<2.5 (Law et al. 
2009; Lehnert et al. 2009; Forster Schreiber et al. 2009). 
The same procedure used for the SDSS galaxies was applied 
to these galaxies. Metallicity is estimated either from R23 
or from [NII]A6584/Hq:, depending on which lines are avail- 
able. AGN are removed using the BPT diagram (Kauffmann 
et al. 2003a) or, when [OIII]A5007 and H/3 are not available, 
by imposing log([NII]A6584/Ha)<-0.3. The [NII]A6584 hne, 
which is usually much fainter than Ha, is not detected in 
several galaxies, but removing these galaxies from the sam- 
ple would bias it towards high metallicities. For these ob- 
jects we have assumed a value of the intrinsic [NIIJA6584 
flux which is half of the upper limiting flux. When neces- 
sary, the published M* have been converted to a Chabrier 



(2003) IMF. For galaxies without observations of both Ha 
and H/3, dust extinction is estimated from SED fitting, and 
we assume that continuum and the emission lines suffer the 
same extinction. In local starburst lines often suffer of higher 
extinctions (j4v(lines)~2.3^v(SED) according to Calzetti 
et al. 2000). We have checked that the inclusion of this ef- 
fect would have little effect on the final relations and on the 
conclusions of this paper. 

Erb et al. (2006) have observed a large sample of 91 
galaxies at z~2.2. Metallicities have been measured only on 
average spectra stacked according to Af*, which has the re- 
sults of mixing galaxies of different SFRs. Despite this prob- 
lem, no systematic differences in metallicity are detected 
with respect to the other galaxies measured individually, 
and the Erb et al. (2006) galaxies are included in the high- 
redshift sample, although without binning them with the 
rest of the galaxies. 



2.3 z=3-4 

A significant sample of 16 galaxies at redshift between 3 and 
4 was observed by Maiolino et al. (2008) and Mannucci et al. 
(2009) for the LSD and AMAZE projects. Pubhshed values 
of stellar masses, line fluxes and metallicities are available for 
these galaxies, which can be compared with lower redshift 
data. The same procedure as at lower redshift was used, with 
the exception that SFR is estimated from H/3 after correc- 
tion for dust extinction, and metallicities are measured by a 
simultaneous fitting of the line ratios involving [OII]A3727, 
n/3 and [OIII]A4958,5007, as described in Maiohno et al. 
(2008). 
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Figure 2. Three projections of the Fundamental MetaUicity Relation among M*, SFR and gas-phase metallicity. Circles without error 
bars are the median values of metallicity of local SDSS galaxies in bin of M* and SFR, color-coded with SFR as shown in the colorbar 
on the right. These galaxies define a tight surface in the 3D space, with dispersion of single galaxies around this surface of ~0.05 dex. 
The black dots show a second-order fit to these SDSS data, extrapolated toward higher SFR. Square dots with error bars are the median 
values of high redshift galaxies, as explained in the text. Labels show the corresponding redshifts. The projection in the lower-left panel 
emphasizes that most of the high-redshift data, except the point at z=3.3, are found on the same surface defined by low-redshift data. 
The projection in the lower-right panel corresponds to the mass-metallicity relation, as in Fig. 1, showing that the origin of the observed 
evolution in metallicity up to z=2.5 is due to the progressively increasing SFR. 



3 THE MASS-METALLICITY RELATION AS A 
FUNCTION OF SFR 

The grey-shaded area in the left panel of Fig. 1 shows the 
mass-metallicity relation for our sample of SDSS galaxies. 
Despite the differences in the selection of the sample and 
in the measure of metallicity, our results are very similar 
to what has been found by Tremonti et al. (2004). The 



metallicity dispersion of our sample, ~0.08 dex, is some- 
what smaller to what have been found by these authors, 
~0.10 dex, possibly due to different sample selections and 
metallicity calibration. The 4th-order polynomial fit to the 
median mass-metallicity relation is: 



12 -I- log{0/H) = 



8.96 -f 0.31m 
-0.017m 



0.23m^ 



'■' + 0.046m'* 



(1) 
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where 7n=log(M*)-10 in solar units. 

We have computed the median metaUicity of SDSS 
galaxies for different values of SFR. Median have been com- 
puted in bins of mass and SFR of 0.15 dex width in both 
quantities. On average, each bin contains 760 galaxies, and 
only bins containing more than 50 galaxies are considered. 
The left panel of Fig. 1 also shows these median metallicities 
as a function of M* . It is evident that a systematic segrega- 
tion in SFR is present in the data. While galaxies with high 
M* (log(Mt)>10.9) show no correlation between metallicity 
and SFR, at low Af* more active galaxies also show lower 
metallicity. The same systematic dependence of metallicity 
on SFR can be seen in the right panel of Fig. 1, where metal- 
licity is plotted as a function of SFR for different values of 
mass. Galaxies with high SFRs show a sharp dependence 
of metallicity on SFR, while less active galaxies show a less 
pronounced dependence. 

A hint of this effect was already noted by Ellison et al. 
(2008), but the different sample selection and the large bins 
in SFR reduced the observed dependence on SFR to a small 
correction of ~0.05 dex with respect to the value determined 
by the mass-metallicity relation. Also Rupke et al. (2008) 
presented evidences for lower metallicities in local galaxies 
with high SFRs, although with a very large scatter. 
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Figure 3. Metallicity dispersion of single SDSS galaxies around 
the FMR. This histogram shows the differences from the median 
computed in bins of 0.05 dex in M^ and SFR. The red line is a 
gaussian with (t=0.053 dex. Positive differences mean metallicities 
of single galaxies larger than the median. 



4 THE FUNDAMENTAL METALLICITY 
RELATION 

The dependence of metallicity on AI^ and SFR can be 
better visualized in a 3D space with these three coordinates, 
as shown in Figure 2. SDSS galaxies appear to define a 
tight surface in the space, the Fundamental Metallicity 
Relation, with metallicity well defined by the values of 
M* and SFR. All the data on this FMR are shown in table 1. 

The introduction of the FMR results in a significant re- 
duction of residual metallicity scatter with respect to the 
simple mass-metallicity relation. The dispersion of individ- 
ual SDSS galaxies around the FMR, shown in Fig. 3, com- 
puted in bins of 0.05 dex in M* and SFR, is ~0.06 dex when 
computed across the full FMR and reduces to ~0.05 dex i.e, 
about 12%, in the central part of the relation where most 
of the galaxies are found. This means that about half of 
the total scatter of the mass-metallicity relation (0.08 dex) 
is due to the systematic effect with SFR, while about half 
is due to intrinsic differences between galaxies and/or un- 
certainties on the measurements. The reduction in scatter 
becomes even more significant when considering that most 
of the galaxies in the sample cover a small range in SFR, 
with 64% of the galaxies (±lcr) is contained inside 0.8 dex. 
Galaxies with very low SFRs are not selected due to the 
high SNR threshold on Ha, while high SFR galaxies are 
rare in the local universe. As a consequence the scatter on 
the full sample is dominated by galaxies having a small cor- 
rection due to SFR. In contrast, considering only galaxies 
with high SFRs, the scatter is reduced by a large factor. For 
example, for log(SFR)>0.5, the scatter around the average 
mass-metallicity relation is about 40%, while it is a factor 
of ~3 lower around the FMR. 

The final scatter is consistent with the intrinsic uncer- 
tainties in the measure of metallicity (~0.03 dex for the 



calibration, to be added to the uncertainties in the line ra- 
tios), on mass (estimated to be 0.09 dex by Tremonti et al. 
2004), and on the SFR, which are dominated by the uncer- 
tainties on dust extinction. Nevertheless, the scatter around 
FMR tends to reduce when the minimum redshift Zmin of 
the galaxy sample is increased. This means that part of the 
residual scatter is probably due to the different apertures 
used to measure mass, based on a total magnitude, and 
metallicity and SFR, derived for the central 3". In partic- 
ular, the effect of metallicity gradients are expected to be- 
come less important at larger redshifts, and for this reason 
the increase of Zmin is able to reduce the scatter. 

A close inspection of fig. 3 reveals the presence of an 
extended wing toward lower metallicities in the distribution 
of scatter. This extension contains ~3% of the objects. 
Most of these galaxies have low M*and high SSFR, and 
could be objects in special conditions. For example, they 
could be interacting galaxies, which will be discussed in 
sec. 7. 

We have fit the median values of metallicity of the SDSS 
galaxies in table 1 with a second-order polynomial in A/* and 
SFR, obtaining: 



12 + log{0/H) 



8.90 + 0.37m -0.14s 
-h0.12ms-0.054s2 



0.19m^ 



(2) 



where ?n=log(M,t)-10 and s=log(SFR) in solar units. The 
residual scatter of median metallicities around this fit is 
0.001 dex. Such a fit, shown in fig. 2, provides a clear 
representation of the dependence of metallicity both on Ah 
and SFR, and allows to compare the local FMR with high 
redshift galaxies. 

The shape of the FMR surface depends on a number 
of factors, such as the selection of the galaxy sample and 
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Table 1. The Metallicity Fundamental Relation for SDSS galaxies selected as in sec. 2.1. For each value of M^, (reported in the first 
line) and SFR (first column) we list the median value of metallicity, the la dispersion around this value, and the number of galaxies in 
each bin. 
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the way metallicity, M* and SFR are measured. We have 
done a number of checks to test whether the result depends 
critically on either of our assumptions. First, we changed 
the thresholds in SNR and redshift used to select the 
galaxy sample, and checked that the results do not change 
systematically with these thresholds. Second, we have 
studied the effect of considering only metallicities derived 



either from R23 or from [NII]A6584/Ha. As discussed in 
sec. 2.1, systematic differences are found but are limited to 
the level of 0.05 dex. Apart from this metallicity offset, the 
shape of the FMR does not change by more than 0.05 dex 
at any point. In particular, there is no large, monotonic 
dependence of the difference with SFR or SSFR. This is 
interesting because there could be systematic effects related 
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either to the density of the HII regions or to the ionization 
parameter U, which could depend on the SSFR (e.g., 
Shapley et al. 2005; Hainhne et al. 2009) and introduce 
spurious behavior. R23 and [NII]A6584/Ha show opposite 
dependencies with U (see, for example, Liu et al. 2008 and 
Nagao et al. 2006). For this reason the absence of systematic 
differences between these two line ratios with SSFR is an 
indication that the dependence of measured metallicity on 
density or U, if present, is small, in agreement with the 
findings of Liu et al. (2008) and Brinchmann et al. (2008). 

There are two points that could affect the shape of the 
FMR. First, SFR is estimated from Hq luminosity corrected 
for extinction using the Balmer decrement. Several authors 
(Kennicutt 1998; Moustakas et al. 2006) have shown that 
this is a reliable SFR indicators over a large range of galaxy 
properties. Others (Chariot & Longhetti 2001; Brinchmann 
et al. 2004) have discussed that systematic effects with mass 
and metallicity could be present. As we have split galaxies in 
bins of mass and SFR, and inside each bin metallicity spans 
a small range, any systematic effect on SFR does not hamper 
the existence of the FMR but could change its shape. Sec- 
ond, we are using SFRs and metallicities that apply only to 
the central 3"of the galaxies, corresponding to 4-11 kpc pro- 
jected angular size given our redshift range. These quantities 
are compared with total mass derived from integrated pho- 
tometry ^ (Kauffmann et al. 2003b; Salim et al. 2007). fn the 
SDSS sample, the fraction of mass contained inside the pro- 
jected fiber aperture (as listed in the MPE/JHU catalog) is 
about 1/3 of the total, with log(total mass) - log(fiber mass) 
= 0.50±0.15, and no systematic dependence on mass. We 
find that total and fiber mass are correlated with SFR and 
metallicity at a similar degree, and we used the total mass as 
this is the quantity usually available in galaxy catalogs. The 
measured metallicity is expected to be a fair representation 
of total metallicity: abundance gradients can be present but, 
at least in local galaxies, they are of modest importance in 
this contest, both because only weak variations with radius 
are usually found (e.g., Magrini et al. 2009) and because 
our apertures sample a significant part of the galaxies. Cor- 
recting fiber SFR into total ones is not straightforward, as 
there is no guarantee that Ha scales radially as luminosity 
or mass. If SFR scales with radius as mass, we expected that 
total SFR are 3 times larger than the fiber ones considered 
here. The use of total SFR would produce a FMR shifted to- 
wards higher SFRs, and its shape could have some changes. 
Nevertheless the small scatter observed in the FMR means 
that the fiber SFR must be well correlated with the total 
one. 

Summarizing, even if the overall shape of the FMR can 
change in different samples of galaxies and depends on sev- 
eral details, the main properties of the FMR are very robust 
and passed all our tests. 



5 THE LOCAL FMR AND HIGH-REDSHIFT 
GALAXIES 

Using the samples described in Sect. 2.2 and 2.3 we compare 
the dependence of metallicity on M* and SFR in galaxies at 
z~0.8, z~1.4, z~2.2 and z~3.3, and compare it with the 
properties of local SDSS galaxies. Galaxies at all redshifts 
follow well defined mass-metallicity relations (see, for ex- 
ample, Mannucci et al. 2009, and references therein). For 
this reason each of these samples, except the one an z~3.3 
that contains 16 objects only, is divided into two equally- 
numerous samples of low- and high-M* objects. Median val- 
ues of Mi,, SFR and metallicities are computed for each of 
these samples. 

Galaxies at intermediate and high redshifts show, on av- 
erage, larger SFR with respect to local SDSS galaxies. This 
is easily explained by selection, because only galaxies with 
significant SFRs are selected and observed spectroscopically. 
Galaxies at z~0.8 have values of M* and SFR which overlap 
with the SDSS sample, and therefore the two samples can 
be directly compared. These galaxies are found to be com- 
pletely consistent with the FMR defined by SDSS galaxies, 
with no evidence for evolution. This is shown in Figs. 2 and 
4. At redshift above 1, some extrapolation towards higher 
SFRs of the fit in eq. 2 in required. All galaxies are within 
0.6 dex from the most active SDSS galaxies, while most mas- 
sive galaxies at z=2.2 require an extrapolation of 1 dex. For 
comparison, SDSS galaxies span two orders-of-magnitude in 
SFR (see Fig. 4). The necessity of an extrapolation intro- 
duces some uncertainty, but we have checked the the result 
does not depend critically on the characteristic of the fit, 
such as the degree of the used polynomial. 

When taking into account the uncertainties, data up to 
z~2.5 are consistent, both in shape and in normalization, 
with the same FMR defined by SDSS, with no evidence 
for evolution. Distant galaxies show larger dispersions than 
the local SDSS galaxies, between 0.2 and 0.3 dex. At least 
part of these relatively larger dispersions are due to the 
large uncertainties in the estimates of both metallicity 
and SFR, but part of it could be intrinsic, related to 
different evolutionary stages of the galaxies. A larger sam- 
ple of well measured galaxies is needed to address this point. 

The existence of a relation between mass, SFR and 
metallicity could be considered a mere consequence of the 
existence of two other well-known relations, the mass- 
metallicity relation, and the mass-SFR relation (see, for 
example, Schiminovich et al. 2007). fn fact, the novelty of 
this relation can be understood by comparing the relative 
redshift evolutions: while both the mass-metallicity and 
the mass-SFR relations are known to evolve significantly 
with redshift (Daddi et al. 2007; Mannucci et al. 2009; 
Forster Schreiber et al. 2009), the FMR remains the same 
up to z=2.5, a period of time spanning 80% of the universe 
lifetime. Under this respect, the FMR seems to be the 
fundamental one, directly related to the mechanisms of 
galaxy formation. 



In the SDSS sample, metallicity changes more with M* 

(~0.5 dex from one extreme to the other at constant SFR, 

see Fig. 1) than with SFR (~0.30 dex at constant mass). 

1 see http://www.mpa-garching.mpg.de/SDSS/DR7/masscomp.html Therefore mass is the main driver of the level of chemical 
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Figure 4. Left: Mctallicity as a function of SFR for galaxies in the three bins of Af* containing high-redshift galaxies. The values of 
log(M*) are shown by the labels on the left. Empty square dots are the median values of metallicity of local SDSS galaxies, with error 
bars showing la dispersions. Lines are the fits to these data. Solid dots are median values for high-redshift galaxies with z<2.5 in the 
same mass bins, with labels showing redshifts. Right: metallicity difference from the FMR for galaxies at different redshifts, color-coded 
in mass as in the left panel. The SDSS galaxies defining the relation are showing at z^O.l with their dispersion around the FMR. All the 
galaxy samples up to z=2.5 are consistent with no evolution of the FMR defined locally. Metallicitics lower by ~0.6 dex are observed at 
z~3.3. 



enrichment of SDSS galaxies. This is related to the fact that 
galaxies with high SFRs, the objects showing the strongest 
dependence of metallicity on SFR (see the right panel of 
fig. 1), are quite rare in the local universe. At high redshifts, 
mainly active galaxies are selected and the dependence of 
metallicity on SFR becomes dominant. 



5.1 Possible effects on data at z>2.5 

Galaxies at z~3.3 show metallicities lower of about 0.6 dex 
with respect to both the FMR defined by the SDSS sample 
and galaxies at 0.5<z<2.5. This is an indication that some 
evolution of the FMR appears at z>2.5, although its size 
can be affected several potential biases of different nature 
that should be taken into careful account. 

First, metallicity at z=3.3 is measured by the oxy- 
gen indicators only, while SDSS galaxies use both R23 and 
[NII]A6584/Hq. As discussed in sec. 2.3, in the local universe 
both indicators give consistent results, with systematic dif- 
ferences limited to 0.05 dex. Also, galaxies at z~0.8 use the 
same metallicity based on R23. Therefore this is not likely 
to be the origin of the evolution of 0.6 dex, although it could 
be responsible for part of the difference. 

Second, systematic evolution with redshift of the pho- 
toionization conditions could be present, for example be- 
cause high-redshift galaxies have larger SFRs. The presence 
of such an effect can be studied by comparing different line 
ratios, and several studies up to z~2.5 indicate that such 
an evolution actually exists (Shapley et al. 2005; Erb et al. 
2006; Brinchmann et al. 2008; Liu et al. 2008; Hainline et al. 
2009) . Nevertheless, this effect is not large enough to move 



galaxies up to z=2.2 off the FMR, even if some of these 
galaxies have SFR larger than most of the galaxies at z=3.3. 
This is in agreement with the results by Liu et al. (2008); 
Brinchmann et al. (2008) who estimate a minor effect of 
this evolution on the measure of metallicity. Therefore evo- 
lution of the conditions could give some contribution to the 
observed evolution but are not likely to be the only reason. 

Third, at z=3.3 the Ha and [NII]A6584 lines fall at 
~2.8/im and are not observed in any galaxy of this group. 
For this reason we cannot use this line ratio to remove AGNs. 
X-ray and mid-IR data on these targets have been analyzed 
in order to exclude dominant AGN (Maiolino et al. 2008), 
but it is possible that some faint AGN is still present among 
these galaxies (Wright et al. 2010). The presence of such 
objects would tend to reduce the measured metallicity. As 
almost all galaxies at z~3 have metallicity below both the 
FMR and the mean values at z~2 (Mannucci et al. 2009), 
this effect can explain the observed difference only if AGNs 
are present in most galaxies. 

Fourth, observations at z~3.3 target the [OIII]A5007 
line, which has a strong dependence on metallicity, with 
more metal-poor regions emitting a brighter [OIIIJA5007 line 
for a given SFR (Maiolino et al. 2008). As only the bright- 
est part of the galaxies are detected in our data, this effect, 
could bias the measured line ratios towards the regions of 
lower metallicities. At lower redshift this effect could be less 
important, both because [OIII]A5007 is usually intrinsically 
fainter than Ha observed at z<2, and because the cosmo- 
logical dimming of the surface brightness, proportional to 
(1 + z)*', is a much more severe problem at higher redshift. 
In our data at z=3.3, metallicity does not seem to increase 
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Figure 5. Left: Residual dispersion of the median values of metallicity of SDSS galaxies as a function of a as defined in eq. 3. The values 
corresponding to the minimum dispersion (a=0.32), to a=0 (/Xa=log(M*)) and and a=l (^a=-log(SSFR)) are shown. Right: metallicity 
as a function of /jo.32, which minimize the residual scatter. Colored lines are local SDSS galaxies, with colors as in the left panel of Fig. 1. 
The black line shows the polynomial best fit. Black dots are high-redshift galaxies, labelled with redshifts. These galaxies, except those 
at z~3.3, appear to follow the same relation define by SDSS galaxies. 



with increasing photometric aperture, although we are lim- 
ited by the low SNR in the external regions of the galaxies. 
This effect could be present but is not likely to produce the 
observed difference of a factor of 3. 

Finally, there could be selection effects resulting in a 
reduction of the average metallicity of observed sample. 
For example, if more metal-rich galaxies also have larger 
amounts of dust, it is possible that our UV-selected galaxies 
at z~3.3 preferentially select low-metallicity objects with 
lower dust column densities. If present, this effect would 
also work at lower redshift, introducing some evolution also 
at z=2.2. Again, such an effect could be present but is 
unlikely to produce a strong differential evolution between 
z=2.2 and z=3.3. 

In conclusion, even if the size of the evolution at z=3.3 
could be affected by several problems, it is unlikely to be 
totally due to observational biases. 



6 A PROJECTION OF THE FMR THAT 

REMOVE SECONDARY DEPENDENCIES 

At a given mass, galaxies with higher SFR have lower metal- 
licities and, therefore, have the metallicity properties of 
lower mass galaxies. As a consequence, we expect that a 
combination of h'h and SFR could be better correlated with 
metallicity. This can be seen in the central panel of fig. 2, 
showing a projection of the FMR which considerably re- 
duces the metallicity scatter. To investigate this point we 
introduce a new quantity fia obtained as linear combination 
of SFR and M* as: 



^Ji^ = log(M.) - Q log(SFR) 



(3) 



where a is a free parameter. For a=0, /io corresponds to 
log(MO, while for a=l, /xi=-log(SSFR). 

The value of a that minimizes the scatter of median 
metallicities of SDSS galaxies around the relation corre- 
sponds to the quantity /ia that is more directly correlated 
with metallicity. Fig. 5 shows the scatter of data in table 1 as 
a function of a. These results show that neither M* (q=0), 
nor, SSFR (a=l) are the quantities producing the smallest 
scatter. In fact, a ~ 0.32 produces a minimum in the disper- 
sion. The resulting diagram is shown in the right panel of 
Fig. 5, where metallicity is plotted against /io.32. The median 
of the resulting distribution cab be fitted by: 

12 + log{0/H) = 8.90 + 0.392:-0.20a;^-0.077a:^ + 0.064x''(4) 

where x = /xo.32 — 10. 

Even if the minimization is computed with SDSS galax- 
ies only, it turns out that high-redshift galaxies up to z=2.5 
follow the same relation between /.10.32 and metallicity as in 
the local universe, and also have the same range of values of 
Po.32. This is the same effect noted in the previous section, 
where high-redshift galaxies have been found to follow the 
extrapolation of the FMR, but with two important changes: 
first, no extrapolation from the SDSS galaxies is now needed, 
because both samples have similar values of /io.32; second, 
it is possible to search for simple physical interpretation of 
/to. 32 in terms of the physical processes in place. 

In practice, metallicity of star-forming galaxies of any 
mass, any SFR and at any redshift up to z=2.5 follow the 
following relation: 



12 + log{0/H) = 



8.90 + 0.472; 
9.07 



ff MO. 32 < 10.2 
if M0.32 > 10.5 



(5) 
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with X = /10.32 — 10. 

6.1 Metallicity and SSFR 

It is interesting to plot metallicity as a function of SSFR, 
as in fig. 6, because several properties of galaxies depend on 
this quantity. In this plot it can be seen that SDSS galax- 
ies of any mass have the same dependence of metallicity 
on SSFR. A threshold SSFR exists, about 10" ^"yr"^ which 
discriminates the abundance effect of the SFR. Above this 
limit, metallicity decreases rapidly with SFR in galaxies on 
any mass. Below this limit, metallicity is constant in massive 
galaxies (log(Aft)>10.5) and slowly decreasing with SFR in 
less massive galaxies. The interpretation of this behavior is 
given in the next section. 

The fraction of galaxies above and below this SSFR 
threshold changes with Af*. Most of the low-mass galax- 
ies, and only a small fraction of high-mass galaxies in our 
SDSS sample are in the "high SSFR" regime, and this is the 
well-known "downsizing" effect (Gavazzi & Scodeggio 1996; 
Cowie et al. 1996). As M* and SFR are independent vari- 
ables of the FMR, downsizing does not shape the relation 
but defines how it is populated, i.e, how many galaxies of a 
given M^: have a certain level of SFR and, as a consequence, 
metallicity. 



T DISCUSSION 

In sec. 4 we have shown that in the local universe a tight 
relation exists between metallicity, stellar mass, and SFR, 
in which metallicity increases with Af* and decreases with 
SFR in a systematic way. In sec. 5 we have shown that the 
same relation, without any evolution, holds up to z=2.5, and 
that the observed evolution of the mass-metallicity relation 
is simply due to the sampling of different parts of this re- 
lation at different redshifts. In sec. 6 we have seen that the 
systematic dependence of metallicity on A/* and SFR at all 
redshifts can be expressed in an easy form by introducing 
the quantity /io.32, linear combination of Mi, and SFR. 

The interpretation of these results must take into ac- 
count several effects. In principle, metallicity is a simple 
quantity as it is dominated by three processes: star forma- 
tion, infall, outflow. If the scaling laws of each of these three 
processes are known, the dependence of metallicity on SFR 
and Mi, can be predicted. In practice, these three processes 
have a very complex dependence of the properties of the 
galaxies, and can introduce scaling relations in many differ- 
ent ways. 

First, it is not known how outflows depend on the prop- 
erties of the galaxies. In many models, star formation pro- 
duces SNe which inject energy, radiation and momentum 
into the interstellar medium, with the result of ejecting part 
of the enriched gas (Veilleux et al. 2005; Spitoni et al. 2010). 
A central Active Galactic Nucleus (AGN) can also provide 
feedback (Somerville et al. 2008). The properties of this 
galactic winds are debated. Dekel & Woo (2003) reproduced 
the mass-metallicity relation with a wind related to the en- 
ergy of SNe, which is proportional to stellar mass. In large 
galaxies with deep potential wells, such a wind is not effec- 
tive in producing an outflow (Tremonti et al. 2004). Murray 
et al. (2005), Dave et al. (2007) and Oppenheimer & Dave 
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Figure 6. Metallicity of SDSS galaxies as a function of SSFR. 
The grey-shaded areas contain 64% and 90% of all SDSS galaxies, 
with the thick central line showing the median relation. The col- 
ored lines show the median metallicities for different values of A/*. 
The lines with black dots show fits to these metallicity distribu- 
tion for four values of log(A^*), 9.4, 9.7, 10, and 10.9. The model 
is described in sec. 7.1 and includes only dilution for infalling gas 
and the SK law. 



(2008) discuss a different scheme, the "momentum-driven 
wind" , in which wind speed increases with galaxy mass while 
its efficiency decreases. 

Second, infalls are expected to influence metallicity in 
two ways. On the one hand, infall of metal poor gas directly 
reduces the observed metallicity by diluting the metal-rich 
gas, as discussed, for example, by Finlator & Dave (2008) 
in the context of their wind scheme. On the other hand, 
infall is expected to produce star-formation activity follow- 
ing the Schmidt-Kennicutt law, and the metals produced by 
new stars are expected to increase metallicity. If merging 
between galaxies, rather than smooth infall from the IGM, 
is the main channel to drive gas into galaxies, the fuel of 
star formation can be significantly enriched, and the dilu- 
tion effect could be absent. 

Third, in some semianalytical models of galaxy forma- 
tion (de Rossi et al. 2007; Mouchine et al. 2008) the be- 
havior of metallicity is dominated by the dependence of 
star-formation efficiency of galaxy mass: less massive galax- 
ies are less evolved and, therefore, show lower metallicities. 
Some other models (Brooks et al. 2007; Dayal et al. 2009; 
Salvaterra et al. 2009) put together downsizing and feed- 
back: metallicity is mainly related to different star formation 
efficiencies in different galaxies, but the efficiency is regu- 
lated by SN feedback. Tassis et al. (2008) agree that low 
star-formation efficiency in low-mass galaxies is the main 
driver of the mass-metallicity relation, but also pointed out 
that mixing of heavy elements in the outer regions of galax- 
ies could help hiding a significant fraction of metals. Also 
Koppen et al. (2007) attribute the different abundances to 
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different levels of productions of metals, but their model in- 
clude systematic variation of IMF, which is proposed to be 
more top-heavy in galaxies with higher SFRs. 

A full exploitation of our results requires a full discus- 
sion of these models, that is beyond the scope of this paper 
and will be subject of a future work. Here we discuss some 
simple implications of our results on the role and the prop- 
erties of infalls and outflows. In the next two subsections we 
will assume two limiting cases for the timescale of chemical 
enrichment as compared to the other relevant timescales. 



detail modelling of this evolution and a comparison of the 
different timescales involved is needed to test if this is a 
viable explanation. 

In this scenario, the dependence of metallicity on SFR is 
due to infall and dominates at high redshifts, where galaxies 
with massive infalls and large SFRs are found. In contrast, in 
the local universe such galaxies are rare, most of the galaxies 
have low level of accretion, and abundances are dominated 
by the dependence on mass, possibly due to outflow. 



7.1 The effect of the infall 

The dependence of metallicity on SFR can be explained by 
the dilution effect of the infalling gas. In a very simple model, 
we assume that galaxies with low SFR have a given gas frac- 
tion fg. When a certain amount Mi„/ of metal-poor gas is 
accreted, galaxies start forming stars at a given SFR de- 
fined by the Schmidt-Kcnnicutt (SK) law on the infalling 
gas. This is an empirical relation between surface densities 
of star formation and surface density of gas: 



-'SFR ' 



(6) 



where n has values around 1.4-1.5 both at low and high 
redshifts (Kennicutt 1998; Bouche et al. 2007; Kennicutt 
2008; Gnedin & Kravtsov 2010; Verley et al. 2010). For SDSS 
galaxies, the spectroscopic aperture is the always same and 
the SK law becomes a relation between masses, obtaining: 



SFR ~ M" 



(7) 



In this model, galaxies are observed during the phase when 
the dilution effect of the infall is predominant over the metal- 
licity enrichment due to new stars, i.e., the observed metal- 
licity [12-|-log(0/H)]oi,s is related to the initial metallicity 
[12+log(0/H)]i„ by the ratio of pre-existing and infalling 
gas: 



[12+log(0/H)]obs = [12+log(0/H)]i„-log 1 + 



Mint 



(8) 



The results of this simple model are presented in fig. 6. For 
simplicity, only 4 values of mass are shown, using the best- 
fitting values of /g=0.3% and Mi^f between 10^'^ and lO'''^ 
M© . This range of values of Mi„f is the same for all masses 
and is determined by the range of observed SSFR. For each 
mass, the metallicity level at low SSFR, [12-|-log(0/H]i„, is a 
free parameter, whose value is possibly determined by other 
effects, such as outflow. Our very simple model is capable of 
reproducing all the main properties of the FMR: metallicity 
reduces with increasing SFR, a threshold of SSFR exists, 
larger metallicity effects are produced in smaller galaxies, 
both above and below the threshold, the slope of the relation 
at high SSFR is exactly as observed. 

For this scenario to work, the timescales of chemical 
enrichment must be longer than the dynamical scales of the 
galaxies, over which the SFR is expected to evolve. In other 
words, galaxies on the FMR are in a transient phase: after 
an infall, galaxies first evolve towards higher SFR and lower 
metallicities. Later, while gas is converted into stars and 
new metals are produced, either galaxies drop out of the 
sample because they have faint Ha, or evolve toward higher 
values of /io.32 and higher metallicities along the FMR. A 



7.2 Properties of outflows 

In many local galaxies, timescales of chemical enrichment 
can be shorter than the other relevant timescales (e.g.. Silk 
1993), and galaxies can be in a quasi steady-state situation, 
in which gas infall, star formation and metal ejection occur 
simultaneously (Bouche et al. 2009). This is the opposite 
situation of what is discussed in the previous section. 

Assuming this quasi steady-state situation, in which 
infall and SFR are slowly evolving with respect to the 
timescale of chemical enrichment, our results can be used 
to derive information on the mechanisms of infall and out- 
flow. Fig. 6 implies that a process exists that depends only 
on SSFR which is effective in reducing metallicity from a 
level that depends on mass. In a steady-state situation, in- 
fall cannot be the only dominant effect. The exponent n of 
the SK relation is larger than 1, and this means that the 
efficiency of star formation increases with gas density, i.e., 
more active galaxies should also be more efficient in con- 
verting metal-poor gas into stars. As a consequence, the SK 
law alone, applied to the infalling gas, would predict the 
opposite to what is actually observed, i.e., metallicities in- 
creasing with SFR at constant mass. Some other effect must 
be present. 

One obvious candidate is the presence of outflows, as 
usually observed in starburst galaxies and discussed by 
many authors. As discussed above, there are several pos- 
sible types of galactic winds, which follow different scaling 
relations with mass and SFR. When the dilution by the in- 
falling gas Minf is considered together with enrichment due 
to the SFR, we can reproduce the dependence of metal- 
licity on /io.32 by introducing an outflow proportional to 
SFR" M*~™, with s and m free parameters. In this case, 
at the first order we obtain: 



12 + log{0/H) ~ log 



SFR 



Minf SFR" m; 



(9) 



and, using the SK law: 



12 + log{0/H) ~ m log(M*) + (1 s) log(SFR) (10) 

where n = 1.5 is the index of the SK law. This is the 
functional form of /Iq, and comparing this equation with 
the best-fitting value of a = 0.32 we obtain m = 1 and 
s — 0.65. In other words, in a steady-state situation the 
outflow must to inversely proportional to mass and must 
increase with SFR°®^ 

We note that, using an index n of the SK relation of 1.5, 
the best-fitting value of a=0.32 corresponds almost exactly 
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to 1 — l/n = 0.33. Using this, the simplest way to combine 
the relevant physical parameters to produce /io.32 is: 

log (^lll^) ~ log(M.) - (1 - l/n)log(SFR) ~ ^0.32 (11) 

In this form, metallicity increases with M*, an efTect easily 
attributable either to downsizing or to outflow. In contrast, 
the dependence of metallicity on Mi„f/SFR is not obvious 
and require complete modelling. 



7.3 Merging and smooth accretion 

The small scatter of SDSS galaxies around the FMR derived 
in sec. 4 can be used to constrain the characteristics of gas 
accretion. For this infall/outflow scenario to work and pro- 
duce a very small scatter round the FMR, two conditions 
are simultaneously required: (1) star formation is always 
associated to the same level of metallicity dilution due to 
infall of metal-poor gas; (2) there is a relation between the 
amount of infalling and outflowing gas and the level of star 
formation. These conditions for the existence of the FMR 
fits into the smooth accretion models proposed by several 
groups (Bournaud & Elmegreen 2009; Dekel et al. 2009), 
where continuos infall of pristine gas is the main driver of 
the grow of galaxies. In this case, metal-poor gas is contin- 
uously accreted by galaxies and converted in stars, and a 
long-lasting equilibrium between gas accretion, star forma- 
tion, and metal ejection is expected to established. 

In contrast, larger scatter around the FMR are expected 
in case of merging, for two reasons. First, in interacting sys- 
tems, part of the gas producing the starburst could be metal- 
poor material due to interaction-induced infall (Rupke et al. 
2008, 2010), but part is expected to be metal-rich material 
already present inside the interacting galaxies. In this case 
the initial dilution is not present, and higher metallicities are 
expected. This is true, in particular, for the smaller mem- 
ber of a merger between galaxies with very different masses. 
The secondary, smaller mass, galaxy is expected to show 
higher metallicity because its star formation activity is ex- 
pected to be fuelled by gas coming from the other, larger, 
more metal-rich, galaxy. This is exactly what is observed by 
Michel-Dansac et al. (2008): higher metallicities, above the 
mass-metallicity relation, are present in the secondary mem- 
ber of minor mergers. Second, the level of SFR is related to 
the properties of the merging galaxies, and is expected to 
vary significantly during the different stages of the interac- 
tions. As a consequence, a large range on SFR is expected 
for a given level of metallicity during the merging history 
of the systems. Both effects are expected to produce larger 
spreads in merging galaxies, and this is what is actually ob- 
served (Kewley et al. 2006; Rupke et al. 2008; Michel-Dansac 
et al. 2008; Peeples et al. 2009). As discussed in sec 4, the 
presence of interacting galaxies in the SDSS could be at the 
origin of the extended wing in the distribution of difference 
with the FMR, and this point will be investigated in a future 
paper. 

The situation at intermediate redshift, up to z=2.5, is 
less clear. Galaxies show larger dispersion in metallicity, but 
this can be explained by larger uncertainties in the measured 
values of metallicity, mass and SFR. With the present data 
sample it is not possible, therefore, to study whether smooth 
accretion is dominant up to z=2.5 as in the local universe or 



if merging has a larger impact on dispersion. Nevertheless, 
the absence of any evolution of the FMR up to these redshift 
support a single physical process of accretion in all these 
galaxies. At even higher redshift, z~3.3, when galaxies show 
large scatter and different average metallicities, it is likely 
that new physical effects become important. 

The uncertainties on SFR can be a critical point. Under 
this respect, the Herschel satellite is expected to improve 
significantly the accuracy of the estimate of total SFRs in 
a short time. With its contribution it will be possible to 
reduce the scatter, especially at high redshifts, study the 
presence of an intrinsic dispersion of properties among the 
galaxies, and obtain a much better characterization of the 
FMR relation at high redshifts. 



7.4 Origin of the evolution at z>2.5 

The interplay between gas accretion, star formation, and gas 
outflow seems to be the same at any redshift up to z=2.5, as 
all these galaxies follow the same FMR. What is important 
in this context is the ratio between the different rates and 
the various timescale involved: gas infall, dynamical times, 
star formation, stellar evolution, supernova explosion, chem- 
ical mixing, outflow. Apparently the relative importance of 
these processes does not evolve up to z=2.5. It is possible 
that at higher redshifts this constant balance does not apply 
any more and lower metallicities are observed while galaxies 
evolve towards the FMR. This could be related to an increas- 
ing importance of merging as a way to drive cold gas into 
the galaxies, at least in the most luminous objects that are 
preferentially selected as LBG. This point can be addressed 
by studying the morphology and the mass-SFR relation in 
these objects. 

It should be noted that the metallicity of galaxies at 
z=3.3 can also be reproduced by the same model in sec. 7.1, 
in which the infall dilution is dominant. Local galaxies are 
reproduced by Mi„/ of the order of 10® — 10^ M©, while gas 
masses of the order of 10^ — 10^° M© are needed for galaxies 
at z=3.3. This is in very good agreement with what obtained 
by Mannucci et al. (2009) and Cresci (2010). Also, galaxies 
at these high redshifts could be preferentially detected dur- 
ing the flrst stages of the starburst, when the dilution effect 
is maximum. This is possible if the starburst has a peak 
on short timescales, shorter than the timescales of metal 
production and chemical mixing, and declines afterwards. 
Galaxies could also become more dust-rich during the later 
phases of the starburst, and drop-out from the UV-selected 
samples. 



8 CONCLUSIONS 

We have studied the dependence of gas-phase metallicity 
12-|-log(0/H) on stellar mass M* and SFR on a few sam- 
ples of galaxies from z=0 to z=3.3. In the local universe, 
we find that metallicity is tightly related to both M* and 
SFR (fig. 1), and this defines the Fundamental Metallic- 
ity Relation (fig. 2). The residual metallicity dispersion of 
local SDSS galaxies around this FMR is about 0.05 dex 
(fig. 3), i.e, about 12%. The well-known mass-metallicity re- 
lation, together with the luminosity-metallicity and velocity- 
metallicity relations, is one particular projection of this rela- 
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tion into one plane, and neglecting the dependence of metal- 
licity on SFR results in doubling the observed dispersion. 

When high redshift galaxies are compared to the FMR 
defined locally, we find no evolution up to z=2.5, i.e., high- 
redshift galaxies follow the same FMR defined by SDSS 
galaxies even if they have higher SFRs (fig 4). This means 
that the same physical processes are in place in the local uni- 
verse and at high redshifts. The observed evolution of the 
mass-metallicity relation is due to the increase of the aver- 
age SFR with redshift, which results in sampling different 
parts of the same FMR at different redshifts. 

At even higher redshift, z~3.3, evolution of ~0.6 dex 
with respect to the FMR is found, although several observa- 
tional effects and selection biases may affect the size of this 
evolution. This is an indication that different mechanisms 
start to dominate. 

Even if the nature of the FMR is 3-dimensional in M*, 
SFR and 12-|-log(0/H), metaUicity is found to be tightly 
correlated with /uo.32 ~ log(M*) — 0.321og(SFR). Galaxies 
at any redshifts up to z=2.5 follow the same /io.32-metallicity 
relation and have the same range of values of /io.32 (fig. 5). 

Metallicity in galaxies of any mass is found to have the 
same dependence of SSFR, with galaxies above the threshold 
of SSFR=10~^*'yr~^ showing a rapidly decreasing metallic- 
ities with increasing SSFR (fig. 6). 

The interpretation of the existence of the FMR, its 
dependence of SSFR, and the role of /io.32 depend on 
the relevant timescales. If dynamical times are shorter 
than timescales for chemical enrichment, the dependence of 
metallicity on SFR can be easily explained by dilution by 
infall. In this case this effect dominates the metallicity evolu- 
tion of galaxies at high redshift, when galaxies grow because 
of massive accretions of metal poor gas and produce large 
SFR. Also galaxies at z=3.3 can fit into this scheme, with 
large masses of infalling gas. This is in agreement with other 
recent independent results (Mannucci et al. 2009; Cresci 
2010). In the local universe, galaxies with large SFR are 
rare and often associated to merging events, and other ef- 
fects becomes dominant which relate metallicity mainly to 
M*. Outfiow is an possibly, although downsizing could also 
work. If, in contrast, infall and SFR evolve on timescales 
much longer than the chemical enrichment timescale, a sort 
of steady-state situation is created: continuos infall of metal- 
poor gas, which both sustains SFR and dilutes metallicity, 
and outfiow of metal-rich gas in galactic winds. In this case 
the outflows must depend on both mass and SFR. 

The small residual scatter around the FMR in the 
local universe supports the smooth accretion scenario, 
where galaxy grow is dominated by continuos accretion of 
cold gas. Interacting and merging galaxies are expected to 
show larger spread around the FMR, in agreement to what 
is actually observed. Galaxies at intermediate and high 
redshifts show larger metallicity dispersions, which could 
be due either to uncertainties in the measurements, or to 
intrinsic dispersion, or both. This effect prevents us to study 
the evolution of residual scatter with redshift. Nevertheless, 
the absence of significant biasses in metallicity or in SFR up 
to z=2.5 points toward the existence of the same physical 
effects and the dominance of smooth accretion even at 
intermediate redshift. 
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